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Abstract 

We develop a 1+1 dimensional hydrodynamical model for central heavy-ion collisions at ultra- 
relativistic energies. Deviations from Bjorken's scaling are taken into account by implementing 
finite-size profiles for the initial energy density. The calculated rapidity distributions of pions, 
kaons and antiprotons in central Au+Au collisions at ^/snn = 200 GeV are compared with ex- 
perimental data of the BRAHMS Collaboration. The sensitivity of the results to the choice of 
the equation of state, the parameters of initial state and the freeze-out conditions is investigated. 
Experimental constraints on the total energy of produced particles are used to reduce the number 
of model parameters. The best fits of experimental data are obtained for soft equations of state 
and Gaussian-like initial profiles of the energy density. It is found that initial energy densities 
required for fitting experimental data decrease with increasing critical temperature of the phase 
transition. 

PACS numbers: 12.38.Mh, 24.10.Nz, 25.75.-q, 25.75.Nq 
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I. INTRODUCTION 



High-energy heavy-ion colhsions provide a unique tool for studying properties of hot 
and dense strongly interacting matter in the laboratory. The theoretical description of such 
collisions is often done within the framework of a hydrodynamic approach. This approach 
opens the possibility to study the sensitivity of collision dynamics and secondary particle 
distributions to the equation of state (EOS) of the produced matter. The two most famous 
realizations of this approach, which differ by the initial conditions, have been proposed by 
La.da„ Q (fun .opp,,,) and B,o..e„ U (partial transparency). In recent decades nranv 
versions of the hydrodynamic model were developed ranging from simplified 1+1 3|, |j, l5|, 



y, 0, 3 and 2+1 dimensional models S 3, 11^ 3 



to more sophisticated 3+1 dimensional 
mention the multi-fluid models {20], I21I 



mode 



24 



15 



Jj o 



16 



the Landau or Bjorken type 



17, 



18 



191 ]. One should also 



25l . l26l | which consider the whole collision 



process including the nuclear interpenetration stage. Recent theoretical investigations show 
that fluid-dynamical models give a very good description of many observables at the SPS 
and RHIC bombarding energies (see e.g. Ref. 2^]). 

The 2+1 dimensional hydrodynamical models have been successfully applied 
, [ll|, [12, [1^ to describe the px distributions of mesons and their elliptic flow at midra- 
pidity. These models assume a boost-invariant expansion [2] of matter in the longitudinal 
(beam) direction and, therefore, cannot explain experimental data in a broad rapidity region, 
where strong deviations from the scaling regime have been observed. More realistic 3+1 di- 
mensional fluid-dynamical simulations have been already performed for heavy-ion collisions 
at SPS and RHIC energies. But as a rule, the authors of these models do not study the 
sensitivity of the results to the choice of initial and final (freeze-out) stages. On the other 
hand, it is not clear at present, which initial conditions. Landau-like pL|] or Bjorken-like ^, 
are more appropriate for ultrarelativistic collisions. 

Our main goal in this paper is to see how well the fluid-dynamical approach can describe 
the RHIC data on vr, K,p distributions over a broad rapidity interval, reported recently by 



the BRAHMS Collaboration 



29( 1 . Within our approach we explicitly impose a constraint 



on the total energy of the produced particles which follows from these data. For our study we 
apply a simplified version of the hydrodynamical model, dealing only with the longitudinal 
dynamics of the fluid. This approach has as its limiting cases the Landau and Bjorken 
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models. We investigate the sensitivity of the hadron rapidity spectra to the fluid's equation 
of state, to the choice of initial state and freeze-out conditions. Modiflcation of these 
spectra due to the feeding from resonance decays is also analyzed. Special attention is paid 
to possible manifestations of the deconflnement phase transition. In particular, we compare 
the dynamical evolution of the fluid with and without the phase transition. 



II. FORMULATION OF THE MODEL 



A. Dynamical equations 

Below we study the evolution of highly excited, and possibly deconflned, strongly inter- 
acting matter produced in ultrarelativistic heavy-ion collisions. It is assumed that after a 
certain thermalization stage this evolution can be described by the ideal relativistic hydro- 



dynamics. The energy-momentum tensor is written in a standard form 53 1 



T^"^ = {e + P)U^'U'' - Pg'"', (1) 

where e, P and are the rest-frame energy density, pressure and the collective 4-velocity 
of the fluid. 

We consider central collisions of equal nuclei disregarding the effects of transverse col- 
lective expansion. It is convenient to parametrize f/^ in terms of the longitudinal flow 
rapidity Y as If^ = (coshy, 0, sinhF)^. All calculations are performed using the light-cone 
variables |2|, namely, the proper time r and the space-time rapidity r] , deflned as 

r = v/t^3^, r] = tanh"^ (-) = -\n ^ . (2) 

In these coordinates, the equations of relativistic hydrodynamics, d^T^^ = 0, for an ideal 
baryon-free fluid take the following form 

(r^ + tanh(r -v)^)e + {e + P) (tanh(r " ^)^^ + >^ = 0, (3) 

(e + P) (r |- + tanh - ^) >^ + (tanh ~ v)r^ + ^) P = ■ (4) 

To solve Eqs. Q-dl]), one needs to specify the EOS, P = P{e), and the initial proflles e(ro, 77) 
and y(ro,?7) at a time t = tq when the fluid may be considered as thermodynamically 
equilibrated. 



Following Ref. 17|, we choose the initial conditions for a finite-size fiuid, generalizing the 
Bjorken scaling conditions: 

(1^1 -^o)^ 



-e{\v\-vo) 



(5) 



where B(x) = (1 + sgna;)/2 . The particular choice i]q = corresponds to the pure Gaussian 
profile of the energy density. At small a such a profile can be similar to the Landau initial 



condition 



54( 1 ■ On the other hand, when a or t]q tends to infinity, one gets the limiting case 



of the Bjorken scaling solution. Below we adopt the value tq = Ifm/ c. 



The numerical solution of Eqs. (l3l)- (HD i s obtained by using the relativistic version 3l|] of 



the fiux-corrected transport algorithm 



32|. 



B. Equation of state 

As well known, a deconfinement phase transition is predicted by quantum chromodynam- 
ics (QCD). This phase transition is implemented through a bag-like EOS in the parametriza- 



tion suggested in Ref. [13|]. This EOS consists of three parts, denoted below by indices H, M 
and Q corresponding, respectively, to the hadronic, "mixed" and quark-gluon phases. In 
the case of equilibrated baryon-free matter the pressure P, energy density e and entropy 
density s may be regarded as functions of the temperature only. The hadronic phase con- 
sists of pions, kaons, meson resonances and baryon-antibaryon pairs. It corresponds to the 
domain of low energy densities, e < en, and temperatures, T < Th- The sound velocity. 



Cs = \/ dP/de, is assumed to be constant (c^ = ch) in this phase: 
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(6) 



The mixed phase corresponds to intermediate energy densities, from eh up to eg. The 
following parametrization is used for this region: 



2 \ D rr I ^~ \ TTT^ 



P = cU-{l + ci)BM, T = Th( ^y-^-M (eH<e<eQ). (7) 

Here is the bag constant, determined from the condition of continuity of -P(e) at e = en ■ 
Due to the small sound velocity cm (see Table [T]), both pressure and temperature increase 
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only weakly with e in the mixed phase region. The third, quark-gluon plasma region of the 
EOS corresponds to energy densities above eg : 

P = cle-{l + cl)BQ,T = TJ^^)^ (e>eQ). (8) 

Here Bq is the bag constant in the deconfined phase. The corresponding formulae for the 
entropy density are obtained from the thermodynamic relation s = (e + P)/T. We use the 
sound velocities cjj, c\,j, Cq close to those used in Refs. 



TABLE I: Parameters of EOSs with the deconfinement phase transition. 















Th 




Bm 


Bq 




(GeV/fm3) 


(GeV/fm3) 








(MeV) 


(MeV) 


(MeV/fm3) 


(MeV/fm3) 


EOS-I 


0.45 


1.65 


0.15 


0.02 


1/3 


165 


169 


-57.4 


344 


EOS-II 


0.79 


2.90 


0.15 


0.02 


1/3 


190 


195 


-101 


605 



The parameters and Tq define the boundaries of a mixed phase region separating 
the hadronic and quark-gluon phases. The critical temperature as defined by lattice 
calculations should lie between Th and Tq , i.e. T^ — (Th+Tq)/2 . Earlier lattice calculations 
(see e.g. Ref. 33|]) predicted the values Tc = (170 ± 10) MeV for the baryon-free two-flavor 
QCD matter. However, a noticeably larger value Tc = (192 ± 11) MeV was reported recently 



in Ref. 



3J] . To probe sensitivity to the actual position of the phase transition, we consider 



two EOSs with different T^ and Tq (see Table H]). The EOS-I corresponds to T^ = 165 MeV 
and the parameters en, eg used in the parametrization LH12 of Ref. In the EOS-II 

we choose Th = 190 MeV and scale en, eg to get the same values of e/T"^ as a function 
of T /Th 3]. Finally, the parameters Bq, Tq are found from the continuity conditions for P 
and T . 

Unless stated otherwise, these EOSs are used in the calculations presented in this paper. 
For comparison, we have performed also calculations with several purely hadronic EOSs. In 
this case we extend Eq. ([6]) to energy densities e > en with the same en^Tn as in Table [H 
but choosing different c^h — 0-15 and 1/3. In Fig. [1] we compare the EOS-I and EOS-II 
as well as two purely hadronic EOSs with constant sound velocities Cg = ch ■ One can see 
that the mixed phase region in the EOS-II occupies larger interval of energy densities, i.e. 
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FIG. 1: Comparison of different EOSs used in this paper. The solid (dashed-dotted) hne is 
calculated using Eqs. ©-([H]) with parameters given by the set EOS-I (II) in Table [H The dashed 
and dotted lines correspond to the hadronic EOSs with constant = 0.15 and 1/3, respectively. 

this EOS has a larger latent heat, eg — en, as compared to the EOS-I. By this reason, the 
life-times of the mixed phase will be longer for the EOS-II, assuming the same initial state. 



C. Total energy and entropy 



Using the equations of fluid dynamics one can show that the total energy and entropy of 
the fluid can be expressed as 3] 

+ 00 

E = j da^.T^'^ = S±To j drj [e cosh Y cosh {Y - rj) + P sinh Y sinh {Y - r/)] , (9) 

— oo 

+ 00 

S = j da^sU^" = S^To Jdris cosh {Y - r]) , (10) 
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where S± is the transverse area of the fluid. The right-hand sides of Eqs. (!9|)- (fT0|) give the 
values of the energy and entropy at r = tq . Equations Qj and ffTOl) can be considered as 
sum rules for the total energy and entropy of the produced particles. 

Below we use Eq. to constrain possible values of the parameters characterizing the 
initial state. This is possible since the total energy of produced particles can be estimated 
from experimental data. Indeed, the value of the total energy loss, AE = 73 ± 6 GeV per 
nucleon, has been obtained from the the net baryon rapidity distribution in most central 
Au+Au collisions 28|]. This gives the estimate of the total energy of secondaries in the 
considered reaction: 

E = iVp,rtAE~26.1TeV, (11) 

where A^part — 357 is the mean number of participating nucleons. Substituting the 
parametrization ([5]) into Eq. (Q and taking the value of E from Eq. (ITT]) , one gets the 
relation between the parameters eo,riQ,a . 

We have checked that our numerical code conserves the total energy E and entropy S at 
any hypersurface cr^ lying above the initial hyperbola r = tq , on the level better than 1% 
up to very long times, r ~ 10^ fm/ c. 

D. Particle spectra at freeze out 

The momentum spectra of secondary hadrons are calculated by applying the standard 



Cooper-Frye formula 35|], assuming that particles are emitted without further rescatter- 
ings from the elements da^ of the freeze-out hypersurface r = Tp{ri) . Then, the invariant 
momentum distribution for each particle species is given by the expression 

where p'^ is the 4-momentum of the particle, y and Pt are, respectively, its longitudinal 
rapidity and transverse momentum, g denotes the particle's statistical weight. The subscript 
F in the collective 4-velocity f/^ , temperature T and chemical potential /i implies that these 



56( 1 . The plus or minus sign in the right- 



quantities are taken on the freeze-out hypersurface 
hand side of Eq. (fT2l) correspond to fermions or bosons, respectively. 

As has been already stated, the effects of transverse expansion are disregarded in our 
approach. Due to this reason, we cannot describe realistically the pt spectra of pro- 



duced hadrons, and analyze below only the rapidity spectra. For a cylindrical fireball with 
transverse cross section S± expanding only in the longitudinal direction, one can write 
da^ = S±{dz, 0, dt)'^ . Using Eq. ([2]) one arrives at the following relation 

d(7i_,p^ = S±mT {TFiv) cosh(?/ - r/) - Tp{i]) sinh(?/ - dr] . (13) 

Here rriT is the particle's transverse mass defined as = y^mM-pf^, where m is the 
corresponding vacuum mass. In the same approximation one can also write the expression 

pMf = cosh {y -Ypir])) , (14) 

where 1f('7) = ^{'^f{v)jV)- explicit expression for particle spectra at freeze-out is 
obtained after substituting ( !T3l) -(fT4 l) into Eq. f fT2l) and integrating over t] from — oo to +oo . 
Note that Bjorken's model corresponds to Yp = rj and Tp,Tp independent of rj. As can 
be seen from Eqs. (|T2l)-( |T^ . the rapidity distributions of all particles should be fiat in this 
case. 

We adopt a freeze-out criterion, assuming that a given fluid element decouples from 
the rest of the fluid when its temperature decreases below a certain value Tp . For finite- 
size initial conditions, T{TQ,r]) ^ at |?7| —>■ oo, so that the fluid elements at large \ri\ have 
temperatures below Tp from the very beginning, i.e. at r = tq. We treat these elements 
as decoupled instantaneously {rp = tq) and use in Eq. (|T2l) the initial values of Y and T 
instead of Yp and Tp . Direct calculation shows, that such elements contribute only little to 
the tails of the rapidity distributions. The value of Tp is considered as an adjustable model 
parameter which is found from the best fit to experimental data. 



E. Feeding from resonance decays 

In calculating particle spectra one should take into account not only directly produced 
particles but also feeding from resonance decays. Below we assume that the freeze-out tem- 
peratures for directly produced particles and corresponding resonances are the same. One of 
the most important contributions to the pion yield is given by p(770)-mesons. The spectrum 



of TT^-mesons originating from these decays is calculated by using the expression [361] 



d^Np^^+ 1 f dmRw{mR) /" 3 d^Nn fppR tur. 



2inn 



8 



where the first integration corresponds to averaging over the mass spectrum of p-mesons, 
Pr and p are, respectively, the 4-momenta of the p-resonance and of the secondary pion. 
The normahzation coefficient in Eq. (1151) takes into account that the number of vr^-mesons 
produced in p-decays equals 2/3 of the total multiplicity of p-mesons. The freeze-out mo- 
mentum spectrum of p-mesons, d^Nji/ d^pn , is calculated using Eqs. (IT^ - (in|) with m = mn, 
g = gp = 9. We use the parametrization of the p-meson mass distribution, wirnR) , sugges- 



ted in Ref. 



36|. 



The feeding of the pion yields from other meson and baryon resonances {R = t], 
uj, K*, A . . .) is obtained in the zero-width approximation, assuming that the contribution of 
the resonance R is proportional to its equilibrium density nji(Tp), multiplied by a factor d^, 
the average number of tt"*" mesons produced in this resonance decay {dp = 2/3, drj = 0.65 . . .). 
The details of riR, and dn, calculations can be found in Ref. [s?]. We have checked for several 
resonances with two-body decays (e.g. for R = K*) that such a procedure yields a very 
good accuracy. As a result, we get the following formula for the total resonance contribution 
to the spectrum of tt^ mesons: 

Ed^NR^T,+ d^Np^T,+ 

where the enhancement factor a is defined as follows 

^dpnp{Tp) ^ ' 

We include meson (baryon and antibaryon) resonances with masses up to 1.3 (1.65) GeV 
and widths F < 150 MeV. The statistical weights, masses and branching ratios of these 
resonances are taken from Ref. 38|] . The factor a decreases gradually with decreasing 
freeze-out temperature: a = 2.8, 2.4, 2.3 for Tp = 165, 130, 100 MeV, respectively. 

When calculating the kaon spectra we explicitly include feeding from decays of i^'*(892) 
(in the zero-width approximation). Higher resonances [R = (p, Ki . . .) are taken into account 
by applying the same procedure as for pions. In this case the enhancement factor changes 
from 1.5 to 1.2 when Tp goes from 165 to 100 MeV. 
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III. RESULTS 



A. Best fits of rapidity spectra 

Below we show the results for rapidity distributions of vr- and A'-mesons as well as 
antiprotons produced in central Au+Au collisions at ^/snn = 200 GeV. In all calculations 
we use the fireball radius R = 6.5 fm and Si = ttE? ~ 133 fm^. The results are compared 



with data of the BRAHMS Collaboration 



28 



29| for 5% most central collisions. 



We have considered different profiles of the initial energy density, ranging from the 
Gaussian-like (r/o = 0) to the table-like (cr = 0). We found that in the case of EOS-1 it is not 
possible to reproduce the BRAHMS data on the pion and kaon rapidity spectra in Au+Au 
collisions by choosing either too small (eo ^ 5GeV/fm^) or too large (eo ^ 15GeV/fm^) ini- 
tial energy densities. For such eg values the pion and kaon yields can not be reproduced with 
any Tp . It is also found that the quality of fits is noticeably reduced for initial energy density 
profiles with sharp edges, corresponding to a < 1 . As follows from the constraint ( fTTl) . such 
profiles should have either very large eo or a wide plateau — ?7o < t] < rjQ. This would lead to 
more flat rapidity distributions of pions and kaons as compared to the BRAHMS data. 

A few parameter sets which give good fits with the EOS-I are listed in Table [Tll AH 
three sets from Table [TT] give very similar rapidity distributions for both pions and kaons. 
In these calculations we choose various eg and a and determine tjq from the total energy 
constraint (fTTl) . It is interesting that the initial states A-C have approximately the same 

TABLE II: Parameters of the initial states which give the best fits of the pion, kaon and antiproton 
rapidity spectra observed in central Au+Au collisions at ^/snn = 200 GeV. All sets correspond 
to the EOS-I. To denotes the maximum temperature at r = tq . Ei and are total energies of 
produced particles within the rapidity intervals |y| < 1 and \y\ < 3, respectively. 



set 


eo (GeV/fm3) 


0" 


Vo 


To (MeV) 


El (TeV) 


E3 (TeV) 


E/S (GeV) 


A 


10 


1.74 





279 


1.53 


9.25 


0.89 


B 


9 


1.50 


0.62 


271 


1.54 


9.59 


0.86 


C 


8 


1.30 


1.14 


263 


1.49 


9.55 


0.86 



total entropy 5* ~ 3 x 10^ . This, in fact, should follow from the correct description of 
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total pion and kaon multiplicities 571] . As one can see from the last column of Table |TTl 
the corresponding E / S-iaXios fall into a narrow interval 0.86 — 0.89 GeV. This observation 
is similar to the result of Ref. [391] that the observed ratio of the rest frame energy to the 
multiplicity of produced hadrons is constant as a function of the bombarding energy. 

To check the sensitivity to the parameters of the phase transition, we also calculate the 
pion and kaon rapidity distributions for the EOS-II. It is found that with the same initial 
energy profiles as for the EOS-1 it is not possible to reproduce the observed spectra at any 
freeze-out temperature. In particular, the predicted kaon yield is strongly overestimated |58|] 
at 100 MeV < Tp < Th = 190 MeV. Nevertheless, the BRAHMS data can be well repro- 
duced with the EOS-II too when taking smaller initial energy densities as compared with 
the EOS-I. Fits of approximately same quality are obtained for ~ 5GeV/fm^. As be- 
fore, in choosing the initial conditions we apply the constraint ffTTj) for the total energy of 
produced particles. Similarly to the case of the EOS-I, the data are better reproduced for 
initial profiles with small rjQ < 1 . 

Figures [2H3] show the model results for pion and kaon rapidity distributions obtained for 
the EOS-I and EOS-II. These results correspond to Gaussian initial profiles with r]o = Q . 



For both EOSs we choose the parameter to obtain the best fit of the BRAHMS data[59|. 
Although the overall fits are very similar for both EOSs, the rapidity spectra obtained 
with the EOS-II are slightly broader than those with the EOS-I. In the same figures we 
demonstrate sensitivity to the choice of the freeze-out temperature. The best fits of the pion 
spectrum for EOS-I and EOS-II are achieved with Tp — 130 MeV (see Fig. [2]). On the other 
hand, the kaon spectrum can be well reproduced only by assuming that kaons decouple at 
the very beginning of the hadronic stage, i.e. at Tp ~ = 165 (190) MeV for EOS-I (II). 
The contribution of resonance decays turns out to be rather significant, especially in the 
central rapidity region, where it amounts to about 35% (45%) of the total pion (kaon) yield. 

According to Fig. [2], larger yields of secondary pions are predicted for smaller freeze-out 
temperatures. A much weaker sensitivity to Tp is found for kaons (see Fig. [3]). This difference 
can be explained by the large difference between the pion and kaon masses. Indeed, in the 
case of direct pions, a good approximation at Tp > 100 MeV is to replace the transverse 
mass mp in Eqs. (fT2|) - f|T^ by the pion transverse momentum pp. Neglecting the second term 
in the right-hand side of Eq. f|T3|) . one can show that the rapidity distribution of pions at 



y = is proportional to = Tp(ri) cosh rj-Tp/ cosh^ ^^(^7) integrated over all rj . For a rough 
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FIG. 2: Rapidity distribution of 7r+-mesons in central Au+Au collisions at ^/snn = 200 GeV. Left 
panel shows results of hydro dynamical calculations for the EOS-I and initial conditions jS]) with 
the parameters eo = lOGeV/fm^, 770 = 0, o" = 1.74 (set A from Table HI]) . Right panel corresponds 
to the EOS-II and the parameters eo = 5GeV/fm^, tjq = 0, a = 2.02. Solid, dashed and dashed- 
dotted curves are calculated for different values of the freeze-out temperature Tp . The dotted 
lines show contributions of resonance decays in the case Tp = 130 MeV. Experimental data are 
taken from Ref. [29 ]. 
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FIG. 3: Same as Fig. [21 but for rapidity distributions, 
estimate, one can use the Bjorken relations [2I Yp = r], spTp = sqTq , where Si;' is the entropy 
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density at T = Tp . Using Eq. ([6]) one gets Tp cc Sp oc Tp " and therefore, C, oc Tp " . 
This shows that for < 1/3 the pion yield grows with decreasing Tp. Quahtatively, one can 
say that at low enough ch the increase of the spatial volume at freeze-out compensates for 
the decrease of the pion occupation numbers at smaller Tp. This effect is somewhat reduced 
because of a decreasing resonance contribution at smaller temperatures. It is obvious that 
for kaons this effect should be much weaker due to the presence of the activation exponent 
exp {—rriK/Tp) . In fact, a numerical calculation for the same EOS and initial state shows 
that the kaon yield changes nonmonotonically: first it slightly increases when temperature 
goes down but then it starts to decrease at lower Tp . 

To study sensitivity of particle spectra to the presence of the phase transition, we have 
performed calculations with purely hadronic EOSs. In this case we use the same initial 
conditions as before and apply Eq. ([6]) for all stages of the reaction, including high density 
states. Our analysis shows that for soft hadronic EOS with cjj ~ 0.1 — 0.2 it is possible 
to reproduce the observed pion and kaon data with approximately the same fit quality as 
in the calculations with the quark-gluon phase transition. Furthermore, the corresponding 
freeze-out temperatures do not change significantly. However, we could not achieve satis- 
factory fits for the "hard" hadronic EOS with cjj ^ 1/3. This is demonstrated in Fig. H] 
where we compare calculations for two EOSs with and without phase transition. In both 
calculations = 1/3. In the first case with we use Eqs. ©-(jH]) with the same T^, €h and 
eg as for the EOS-I, but choose cjj- = 1/3. The hadronic EOS is obtained by extending 
Eq. ([6]) to all energy densities. We have found that calculations with cjj = 1/3 require much 
higher initial energy densities as compared to the EOS with cjj = 0.15. One can see that 
this hadronic EOS predicts a too wide pion rapidity distribution. The same conclusion is 
valid for kaons. The reason is that the higher pressure gives a stronger push to the matter 
in forward and backward directions. From these findings we conclude that a certain degree 
of softening of the EOS is required to reproduce the pion and kaon rapidity distributions. 

It turns out that our model can also reproduce reasonably well the antiproton rapid- 
ity spectra measured by the BRAHMS Collaboration 28|] . Figure |5] shows the antiproton 
rapidity distributions, calculated for the EOS-I and the parameter set A. In this case we 
explicitly take into account the contribution of the A(1232) vrp decays, ignoring the width 
of A-isobars. Contributions of higher antibaryon resonances are taken into account in a sim- 
ilar way as for pions and kaons. The resonance contribution is about 55% at Tp = 165 MeV. 
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FIG. 4: Rapidity distributions of tt "'"-mesons in central Au+Au collisions at ^/sNN = 200 GeV 
calculated for the initial condition with the parameters eo = 50GeV/fm'^, rjQ = 0, a = 0.96. Solid 
line corresponds to the EOS with the phase transition (see the text). The dashed line is calculated 



the purely hadronic EOS. In both cases = 1/3. Experimental data are taken from Ref. 



291]. 



One should consider these results as an upper bound for the antiproton yield. A more 
realistic model should include effects of the nonzero baryon chemical potential which will 
cert ainly reduce the antibaryon yield. The thermal model analysis of RHIC data, performed 
in Refs. 



40 



411 ]. gives rather low values for the baryon chemical potentials, /if ~ 30MeV, 



at midrapidity. This will suppress the antiproton yield by a factor ~ exp (— /x^/Tp) ~ 0.8 



B. Rapidity distribution of total energy 

We have calculated additionally the rapidity distribution of the total energy of secondary 
particles, dE/dy, in order to check the energy balance in the considered reaction. In this 
calculation we take into account not only direct pions and kaons (charged and neutral), but 
also heavier mesons and BB pairs (the same set of resonances as in the calculation of pion 
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FIG. 5: Rapidity distributions of antiprotons in central Au+Au collisions at ^/sJfN = 200 GeV. 
Shown are results of hydro dynamical calculations for the EOS-I and the parameter set A. All 
results are obtained assuming fip = 0. Experimental data are taken from Ref. 28 1. 



and kaon spectra). The contribution of heavy mesons and BB pairs was calculated in the 
zero-width approximation at the temperature Tp = 165 MeV. By integrating dE/dy, we 
have determined Ei and E3, the total energies of secondaries within the rapidity intervals 
\y\ < 1 and \y\ < 3, respectively. The BRAHMS Collaboration estimated Ei^^ from the ra- 
pidity distributions of charged pions, kaons, protons and antiprotons in most central Au+Au 
collisions at ^/snn = 200 GeV. The values Ei ~ 1.5 TeV, E3 ^ 9TeV have been reported 



in Ref. 421]. From Table HTl one can see that these values are well reproduced by the model. 

Based on the above analysis we conclude that within the hydrodynamical model the 
BRAHMS data can be well described with the EOS-I and EOS-H and the parameters of the 
initial state (tq = Ifm/c) cr ~ 1.5 — 2, i]o ^ 1 ■ The maximal initial energy density, eo, is 
sensitive to the critical temperature of the phase transition. For the EOS-I {T^ ~ 167 MeV) 
we get the estimate eo — 9 ± 1 GeV/fm^ while for the EOS-II (Tc ~ 192 MeV) the required 
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values of eo are lower by about a factor of two. 

These profiles are intermediate between the Landau and Bjorken limits. It is worth noting 
that the observed pion rapidity distribution can be well approximated by the Gaussian with 



the width o"cxp — 2.3 



29 



given by the expression 



According to the Landau model, the width of the distribution is 

2 



2 



8 \/snn 



where m-Ar is the nucleon mass. For = 1/3 this gives crLan — 2.16 , the value often quoted 
in the literature (see e.g. [29]). On the other hand, for = 0.15 (which is preferable 
within our model) the width is only 1.38 i.e. noticeably smaller than observed by the 
BRAHMS Collaboration. This shows that deviations from the simple Landau model are 
rather significant. 



C. Dynamical evolution of matter 
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FIG. 6: Time dependence of temperature as a function of rj calculated for the parameter set A 
(only the forward hemisphere is shown). Left and right panels correspond, respectively, to the 
EOS-I and the hadronic EOS P = c%e with c% = 0.15 . 

Finally, after we have determined the initial conditions which lead to a reasonable de- 
scription of the observed rapidity spectra, we can use the strength of the hydrodynamical 
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model to follow the dynamical evolution of matter. Below we present the results for two 
EOSs, with and without the phase transition, for the parameter set A. Figures [6H7] show 
profiles of the temperature and the collective rapidity at different proper times r. The main 
difference is that in the case of phase transition the model predicts appearance of a flat 
shoulder in T{ri) and local minima in Y{r]) , which are clearly seen at r < lOfm/c. This 
is a consequence of the mixed phase which has a life time At ~ lOfm/c. According to 
Figs. [6H71 the "memory" of the quark-gluon phase is practically washed out at r > 30 fm/ c . 
As one can see from Fig. [TJ at such times deviations from the Bjorken scaling {Y = rf) do 
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rj 7] 
FIG. 7: Same as Fig. [6l but for collective rapidity profiles. 

not exceed 5%. 

Figure [8] shows the matter isotherms in the rj — t plane. One can clearly see that the 
initial stage of the evolution, when matter is in the quark-gluon phase, lasts only for a 
very short time, of about 5 fm/c. The region of the mixed phase is crossed in less than 
10 fm/c. This clearly shows that the slowing down of expansion associated with the "soft 
point" of the EOS plays no role, when the initial state lies much higher in energy density 
than the phase transition region. In this situation the system spends the longest time in the 
hadronic phase and the late evolution is not sensitive to the phase transition. The freeze-out 
at Tp = 130 MeV requires an expansion time of about 60 fm/c at r/ = 0. This is certainly 
a very long time which is apparently in contradiction with experimental findings. Indeed, 
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FIG. 8: Isotherms in the r] — t plane calculated for the parameter set A. Left and right panel 
corresponds to the same EOSs as in Fig. [6]. Shaded region indicates the mixed phase. 



the interferometric measurements [4^ show much shorter times of hadron emission, of the 
order of 10 fm/ c. This discrepancy can not be removed by considering other EOS or initial 
conditions. A considerable reduction of the freeze-out times can be achieved by including 



the effects of transverse expansion and chemical nonequilibrium 



171]. However, this will 



not change essentially the dynamics of the early stage (r < lOfm/c) when expansion is 
predominantly one-dimensional. A more radical solution of the "short time puzzle" could 
be an explosive decomposition of the quark-gluon plasma, proposed in Ref. 4a]. This may 
happen at very early times, right after crossing the critical temperature line, when the 
plasma pressure becomes rather small. We shall consider this possibility in a forthcoming 
publication. 



IV. SUMMARY AND DISCUSSION 



In this paper we have generalized Bjorken's scaling hydrodynamics for finite-size profiles 
of energy density in pseudorapidity space. The hydrodynamical equations were solved nu- 
merically in r — ?7 coordinates starting from the initial time tq = 1 fm/ c until the freeze-out 
stage. The sensitivity of the final particle distributions to the initial conditions, the freeze- 
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out temperature and the EOS has been investigated. A comparison of vr, K, p rapidity 
spectra with the BRAHMS data for central Au+Au colhsions at ^/snn = 200 GeV has been 
made. Best agreement with these data is obtained for initial states with nearly Gaussian 
profiles of the energy density. In choosing the initial conditions we impose the constraint 
on the total energy of produced particles. It is found that the maximum energy density of 
the initial state, eo, is sensitive to the parameters of a possible deconfinement transition. 
The BRAHMS data are well reproduced with eo of about 10 (5) GeV/fm^ for the critical 
temperature Tc ~ 165 (190) MeV. The only unsatisfactory aspect of these calculations is the 
prediction of a very long freeze-out times, ~ 50 fm/ c for pions. 
We would like to comment on several points. 

It is clear that our 1+1 dimensional model can not be valid at late stages of heavy-ion 
collisions, and the transverse flow effects should be included into a more realistic approach. 



12,ll3|, which 



On the other hand, the above-mentioned 2+1 dimensional models j^l, [lo|. 111 , 
assume Bjorken scaling in the beam direction, are apparently not accurate too, even for 
the slice around 1] = 0. Indeed, in contrast to the Bjorken model, our calculations show 
that total entropy in different pseudorapidity intervals does not stay constant during the 
expansion. Due to the pressure gradients along the beam axis and corresponding fluid's 
acceleration, the entropy is transferred from central pseudorapidity bins to the periphery. 
For instance, the entropy in the central bin |?7| < 1 drops by about 15% during the evolution. 
Therefore, only full 3D models can provide a more reliable description. 

It is interesting to note that the viscosity terms, omitted in this paper, should lead to 
the opposite effects, namely to slower cooling and smaller acceleration of the fluid [46|, |47[ . 
Therefore, to describe the observed data, we would need somewhat broader initial energy 
density profiles and accordingly lower eo values. In principle, our simple model can be used 
for a more quantitative study of these effects. 

We have performed calculations with the initial time To = Ifm/c. Of course, one can 
start the hydrodynamical evolution from an earlier time, i.e. assuming smaller tq . In this 
case one should choose accordingly higher initial energy densities. But To cannot be taken too 
small, since at very early times the energy is most likely stored in strong chromofields [48 1. 
The quark-gluon plasma is produced as a result of the decay of these fields (see e.g. Ref. [49 1 
and references therein) . Estimates show that the characteristic decay times are in the range 
0.3 — 1.0 fm/ c. At earlier times the system will contain both the fields as well as produced 
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partons, and the evolution equations will be more complicated, see e.g. Refs. |50l. |51|. 

It is obvious that the Cooper-Frye scenario of the freeze-out process, applied in this 
paper is too simplified. This was demonstrated e.g. in Ref. js^]- One should also have in 
mind that the freeze-out temperatures obtained in our model will be modified by the effects 
of transverse expansion and chemical nonequilibrium. Attempts to achieve a more realistic 
description of the freeze-out stage have been recently made in Refs. [lo|, 3] where a transport 
model was applied to describe evolution of the hadronic phase. In this approach the solution 
of fluid-dynamical equations is used to obtain initial conditions for transport calculations 
at later stages of a heavy-ion collision. We are planning to use a similar approach in the 
future. 
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Units with h = c = 1 are used throughout the paper. 

Within the Landau model [l| cr oc and eg oc 7^ where 7 is the cm. Lorentz-factor of 
colliding nuclei. 

This is achieved by choosing the same £i/T^ (i = H,Q) for these two EOSs. 

Below we assume that the chemical and thermal freeze-out hypersurfaces coincide. In this 

case /ii? = for baryon-free matter. 

According to the Landau model [1], the multiplicity of produced particles is proportional to 
the total entropy. 

We have checked that at fixed Tp and the same initial conditions the freeze-out times predicted 
for the EOS-II are noticeably longer than for the EOS-I. 

We did not try to achieve a perfect fit of these data, bearing in mind that their systematic 
errors are quite big, about 15% in the rapidity region \y\ > 1.3 29|]. 
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